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Abstract. A search for the axioelectric absorption of 5.5-MeV solar axions produced in the p + d — > 
3 He + 7(5.5MeV) reaction was performed with two BGO detectors placed inside a low-background setup. 
A model independent limit on axion-photon and axion-nucleon couplings was obtained: \gAe x Qan\ < 
3.2 x 10 -9 (m,A = 0). Constraints on the axion-electron coupling constant were obtained for axions with 
masses in the (0.1 — 1.0) MeV range: gAe < (1.8 — 9.0) x 10~ 7 . The solar positron flux from A — > e~ + e + 
decay was determined for axions with masses ttia > 2m e . Using the existing experimental data on the 
interplanetary positron flux, a new constraint on the axion-electron coupling constant for axions with 
masses in the (1.2 — 5.4) MeV range was obtained: gAe < (1 — 5) X 10~ 17 . 

PACS. 1 4.80.Mz,29.40.Mc, 26.65. +t 



1 INTRODUCTION 

The axion concept was introduced in a theory by Weinberg 
[I] and Wilczek [2J, who showed that the solution to the 
problem of CP conservation in strong interactions, pro- 
posed earlier by Peccei and Quinn [3J, should lead to the 
existence of a neutral pseudoscalar particle. The original 
WWPQ axion model contained certain strict predictions 
for the coupling constants between an axion and photons 
(gA-y), electrons (gAe), and nucleons ($an), and was soon 
disproved by experiments performed on reactors and ac- 
celerators, and by experiments with artificial radioactive 
sources [J. 

Two classes of new theoretical models of an "invisi- 
ble" axion retained this particle in the form required for 
solving the CP problem of strong interactions and at the 
same time suppressed its interaction with matter. These 
were the models of "hadronic" (or KSVZ) axion [5], [5] and 
DFSZ axion [TJ,[8] . The scale of Peccei-Quinn symmetry 
violation (Ja) in both models is arbitrary and can be ex- 
tended to the Planck mass mp « 10 19 GeV. The axion 
mass in these models is determined by the axion decay 
constant /a' 

m A (cV) « Umr/fA » 6.0 x 10 6 //a(GcV), (1) 

where w, and /„• are, respectively, the mass and decay 
constant of the neutral 7r meson. Since the axion-hadron 
and axion-lepton interaction amplitudes are proportional 
to the axion mass, the interaction between an axion and 
matter is suppressed. 
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The results from present-day experiments are inter- 
preted within these two most popular axion models. The 
main experimental efforts are focused on searching for an 
axion with a mass in the range of 10~ 6 to 10~ 2 cV. This 
range is free of astrophysical and cosmological constraints, 
and relic axions with such mass are considered to be the 
most likely candidates for the particles that form dark 
matter. 

The effective coupling constants gA-y, gAe, and gAN 
are to a great extent model dependent. For example, the 
hadronic axion cannot interact directly with leptons, and 
the constant g Ae exists only because of radiative correc- 
tions. The constants gA-y and gAN can differ by more than 
two orders of magnitude from the values accepted in the 
KSVZ and DFSZ models [5]. 

New possibilities for solving the CP problem are based 
on the concept of the existence of a world of mirror parti- 
cles [TU] and supersymmetry [TT]. These models allow for 
the existence of axions with a mass of about 1 MeV, and 
this existence is precluded by neither laboratory experi- 
ments nor astrophysical data. 

The purpose of this study is to search experimentally 
for solar axions with an energy of 5.5 MeV, produced in 
the p + d — » 3 He + A (5.5 MeV) reaction as a result of 
the axioelectric effect in bismuth atoms. The axion flux is 
proportional to the pp-neutrino flux, which is known with 
a high degree of accuracy |12) . The range of axion masses 
under study was expanded to 5 MeV. 

In previous works we searched for solar axions emitted 
in the 478 keV Ml-transition of 7 Li [13], in the 14.4 keV 
Ml-transition of 57 Fe [1411T5] and axions produced by the 
inverse Primakoff conversion of photons in the electric field 
of the plasma [151117] . 
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The results of laboratory searches for the axion as well 
as the astrophysical and cosmological axion bounds one 
can find in [¥]. 



2 AXION PRODUCTION IN NUCLEAR 
MAGNETIC TRANSITIONS AND THE 
AXIOELECTRIC EFFECT 

If the axion does exist, the Sun should be an intense source 
of axions. They could be produced efficiently on the Sun 
due to photon-axion conversion in an electromagnetic field 
of plasma. Monochromatic axions would be produced in 
magnetic transitions in nuclei whose low-lying levels were 
excited by the high temperature of the Sun. 

The reactions of a main solar cycle could also produce 
axions. The most intensive flux would be expected as a 
result of the formation of the 3 He nucleus: 



P 



5 He + 7(5.49 MeV). 



(2) 



According to the standard solar model (SSM), 99.7% 
of all deuterium is produced as a result of the fusion of 
two protons, p + p — > d + e + + z/ e , while the remaining 
0.3% is due to the p + p + e~ — > d + v e reaction. The pro- 
duced deuteron captures a proton with lifetime r = 6s. 
The expected solar axion flux can thus be expressed in 
terms of the pp-neutrino flux, which is 6.0 x 10 10 cm _2 s _1 
[12] . The proportionality factor between the axion and 
neutrino fluxes is determined by the axion-nucleon cou- 
pling constant gAN, which consists of isoscalar g AN and 
g AN isovector components. The ratio of the probability 
of a nuclear transition with axion production (uja) to the 
probability of a magnetic transition (w 7 ) takes the form 



Ua 



1 



1 



2na 1 + S 2 



9anP + 9 AN 



(Ho - 0.5)/3 + fi 3 - r\ 



2 / x 3 
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where p 7 and pa are, respectively, the photon and axion 
momenta; 8 2 — E/M is the probability ratio for the E and 
M transitions; a w 1/137 is the fine-structure constant; 
Mo = Mp + Mn ~ 0-^8 and [13 — /i p — /i n « 4.71 are, 
respectively, the isoscalar and isovector nuclear magnetic 
momenta; and /3 and 77 are parameters dependent on the 
specific nuclear matrix elements. 

Within the hadronic axion model, the constants g AN 
and g AN can be written in terms of axion mass [5], [2 1|: 



g° AN = -4.03 x 10- 8 (m A /leV), 



3 
9an 



-2.75 x 10" 



- 8 (m A /leV). (5) 

The similar relations for g AN and g AN for the DFSZ axion 
are model dependent to a greater extent but have the same 
order of magnitude. Their numerical values lie in the range 
of (0.3 — 1.5) from the values of these constants for the 
hadronic axion. 

In the p+d —¥ 3 Hc+7 reaction, the Ml-type transition 
corresponds to the capture of a proton with zero orbital 
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Fig. 1. Ratio of emission probabilities for axions and 7 quanta 
(lja/wj) in the p + d — > 3 He + 7 reaction (curve 1, left-hand 
scale); cross section of the axioelectric effect for 5. 5- MeV axions 
on bismuth atoms for gAe = 1 (curve 2, right-hand scale). 



momentum. The probability of proton capture from the 
S state at proton energies below 80 keV was measured in 
[2"2"] ; at a proton energy of ~ 1 keV, it is % = 0.55 (S 2 — 
0.82). The proton capture from the S state corresponds 
to the isovector transition, and the ratio w^/w 7 , which is 
(3) determined by expression (3), therefore depends only on 
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The calculated values of the uja/u 7 ratio as a function 
of the axion mass are shown in FigJTJ The axion flux on 
the Earth's surface is 

$A = fvppiuA/Ui) = 2.7 x l0 10 (g AN ) 2 { PA /p 7 ) 3 = 

= 2Mxl0- 5 m A (pA/p 7 ) 3 (7) 

(4) where <S> V n V — 6.0 x 10 10 cm _2 s _1 is the pp-neutrino flux 



and rriA is the axion mass in eV units. 

To detect 5.5-MeV axions, we chose the reaction of ax- 
ioelectric effect -i|Z|e-)Z + e. The cross section for 
bismuth atoms exceeds that of the Compton conversion of 
an axion (A + e~ — > 7 + e~) by almost two orders of mag- 
nitude; the detection efficiency for the produced electron 
is close to 100%; and, finally, the background level at 5.5 
MeV is much lower than in the range of natural radioac- 
tivity. As a result, the sensitivity to constants gAe and 
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gAN can be high even in an experiment with a relatively 
small target mass. 

In the axioelectric effect (an analog of the photoelec- 
tric effect), an axion disappears and an electron with an 
energy of E e = Ea — Eb , where Eb is the electron binding 
energy, is emitted from the atom. The axioelectric effect 
cross section for K-shell electrons was calculated (on the 
assumption that Ea >• Eb and Z <C 137) in [25] : 



OAe = 2(Zam) 
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The dependence of the cross section on the axion mass 
for the coupling constant gAe = 1 is shown in Fig. 1. 
The cross section depends on the nuclear charge accord- 
ing to the Z 5 law, and it is therefore reasonable to search 
for this process using detectors with large Z. The K-shcll 
electrons make the main contribution to the cross section. 
The contribution from the other electrons was incorpo- 
rated by introducing a factor of 5/4, by analogy with the 
photoelectric effect. 



3 INTERACTION OF AXIONS WITH 
SOLAR MATTER AND AXION DECAYS 

The flux of 5.5 MeV axions on the Earth-s surface is pro- 
portional to the pp-neutrino flux only when the axion life- 
time exceeds the time of flight from the Sun and when the 
flux is not reduced as a result of axion absorption by solar 
matter. The axions produced at the center of the Sun must 
pass through a layer of ~ 7 x 10 35 electrons per cm~ 2 in 
order to reach the Sun's surface. The Compton conversion 
of an axion into a photon imposes an upper limit on the 
sensitivity of Earth-bound experiments to the constant 
gAe- The cross section of this reaction for 5.5-MeV axions 
depends weakly on the axion mass and can be written as 
a cc w g Ae 4: x 10~ 25 cm 2 . For g Ae values below 10~ 6 , the 
axion flux is not substantially suppressed. 

As can be seen in Fig. 1, the maximum cross section of 



the axioelectric effect on bismuth atoms is GAe ~ 9Ae^-^ x 
10 _23 cm 2 . The abundance of heavy (Z > 50) elements in 
the Sun is ~ 10 -9 in relation to hydrogen. If gAe < 10~ 3 , 
the change in the axion flux does not exceed 10%. The 
axion-photon interaction, as determined by the constant 
gAj, leads to the conversion of an axion into a photon in a 
nuclear field. The cross section of the reaction on a proton 
10 _29 cm 2 , and the condition for the 



is Or 



9a^2 



axion emission from the Sun imposes a constraint on gA-y '■ 
gAj < 10 _4 GeV . The constraint for the other elements 
does not grow stricter, due to their low concentration in 
the Sun. 

The axion-nucleon interaction leads to axion absorp- 
tion in a reaction similar to photo-dissociation: A + Z — , 




10- 17 -J 



m , MeV 



Fig. 2. Constraints on gAe and toa from A — ¥ e + + e~ decay. 
Curves 1 and 2 were obtained in accordance with the condi- 
tion that 90% of all solar axions reach the Sun's surface and 
the Earth's surface, respectively. Curve 3 corresponds to the 
observed positron flux entering the Earth's atmosphere. The 
allowed gAe and vha values are below the corresponding curves. 



Z\ + Zi- This reaction is of a threshold nature; for 5.5- 
MeV axions, it might occur only for a few nuclei: 17 0, 13 C, 
and 2 H. It was shown in [33] that axio-dissociation cannot 
substantially reduce the axion flux for gAN < 10 . 

For axions with a mass above 2m e , the main decay 
mode is related to their disintegration into an electron- 
positron pair: A — , e + + e~ . The lifetime of an axion in 
the intrinsic reference system has the form: 



8tt/(<7 



4m 2 ). 



The probability of an axion reaching Earth is 
P{m A ,PA) = exp(-r os /r cm ), 



(9) 



(10) 



where r as = Lm a/ ' cp a is the time of flight in the reference 
system associated with the axion and L = 1.5 X 10 13 cm 
is the distance from the Earth to the Sun. The condition 
r as < 0.1r CTO (in this case, 90% of all axions reach Earth) 
yields the sensitivity limits for the constant gAe in our 
experiment (Fig. 2, curve 2). 

The positron flux with an energy of 5 MeV near the 
surface of the Earth's atmosphere surface is ps 10~ 4 cm~ 2 sr~ 
Using this value and expression (5) for wa/w-p we can 
obtain the upper constraint on the coupling constant gAe 
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(Fig. 2, curve 3). Curve 1 in Fig. 2 is plotted in accordance 
with the condition that 90% of all solar axions emerge 
from the Sun. The range of excluded gAe and uia values 
thus lies between curves 1 and 3. The obtained upper lim- 
its on the constant gAe {gAe < 10~ 16 — 10~ 17 ) for axions 
with mass in the (1.2 — 5.4) McV range are the strongest 
up-to date. 

If the axion mass is less than 2m e , A — > e + + e~ decay 
is impossible, but the axion can decay into two 7 quanta. 
The probability of decay, which depends on the axion- 
photon coupling constant and the axion mass, is given by 
the expression 



TA^~ n = 647r/^ 7 r 



(11) 



The present-day experimental constraint on gA 7 is 10 



GeV , which corresponds to r c . 



KF s for 1-McV ax- 



ions. This means that the axion flux is not reduced prac- 
tically due to the A — >• 2j decay for axions with masses 
below 1 MeV. 



4 EXPERIMENTAL SETUP 
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We used two scintillation detectors (BGOl and BG02), 
manufactured from orthogermanate bismuth B13G4O12, 
to search for the 5.5 MeV axions. The BGO crystal, which 
had a mass of 580 g, was shaped as a prism 50 mm in 
height, with a cross section in the form of a regular hexagon 
inscribed in a circle 50 mm in diameter. The detector sig- 
nal was measured by an PMT-176 photoelectron multi- 
plier, which had an optical contact with a crystal end sur- 
face. 

External 7 activity was suppressed using passive shield 
that consisted of successive layers of lead (50 mm), iron 
(35 mm), copper (10 mm), and bismuth (15 mm Bi20,3). 
The total thickness of the passive shield was w 100 g cm~ 2 . 
The setup was located on the Earth's surface. To suppress 
the cosmic-ray background, an active shielding composed 
of five plastic scintillators 500 x 500 x 120 mm in size was 
used. The total count rate from the active shielding was 
set at a level of 600 s^ 1 , which led to 4% dead time at an 
inhibit pulse width of 70 /is. 

The spectrometric channel of the BGO scintillation de- 
tector included an amplifier with a shaping time of 1 /is 
and a 12-digit ADC. The amplification was selected so 
that the ADC channel scale was 3.25 keV. Standard cali- 
bration sources ( 60 Co and 207 Bi), in combination with the 
natural radioactivity lines of 40 K and the uranium and 
thorium families, were used for energy calibration of the 
detector. The energy dependence of detector resolution a 
can be presented as a/E s=s 6.3% x _E -1 / 2 , where E is in 
MeV. The range of 5.5-McV electrons in the BGO crys- 
tal was w 3 mm, leading to a loss of electron detection 
efficiency near the detector surface. Based on this elec- 
tron range, our estimate of the detection efficiency was 
e = 0.87. 



Fig. 3. The energy spectrum of the BGOl detector measured 
(1) in anticoincidence and (2) in coincidence with the active 
shielding signal. The arrows indicate peaks with energies of 
1.460 MeV, 2.614 MeV, and 5.5 MeV (desired peak). 



5 RESULTS 

The measurements were performed over 65 days in real 
time by series, each of which was 2 h long. The measure- 
ments were divided into series so as to monitor the time 
stability of the BGO detector and the active shielding. 
The energy spectrum of the BGOl detector in the range 
of 1-13 MeV is shown in Fig. 3. The spectrum of the BGOl 
signals that were not accompanied by the active shielding 
signal is designated as 1. 

In the spectrum, we can identify two pronounced peaks 
at 1.460 and 2.614 MeV; these are due to the natural ra- 
dioactivity of the 40 K (located in the PMT's glass housing) 
and of 208 T1 from the 232 Th family (Fig. 3). The positions 
and intensities of these peaks were used for monitoring of 
time stability. 

Bismuth has the largest nuclear charge among the sta- 
ble isotopes (Z = 83), and the cross section of (e+e~ ^pah- 
production upon the interaction of 7 quanta is therefore 
the largest for this element. The annihilation peak at 0.511 
MeV is pronounced in the spectrum. The peak at 2.1 MeV 
is related to the emission of one annihilation 7 quantum 
from the detector upon the detection of 2.614 MeV 7 rays. 
The kink at 7.5 MeV is due to the 7 quanta produced as 
a result of the capture of thermal neutrons by the nuclei 
of the iron, copper, and lead components of the passive 
shield. 
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Fig. 4. Results from fitting the BGOl spectrum in the range of 
(4.5 — 6.5) MeV. Curve 2 corresponds to the detector response 
function {E = 5.49 MeV, a = 0.148 MeV). 



The positions and dispersion of the 1.46 and 2.614 
MeV peaks determined during the measurements were 
used to find the energy scale and energy dependence of the 
BGO resolution. The energy calibration of the spectromet- 
ric channel was found as a linear function: E — Ax N + B, 
where E is the released energy and N is the channel num- 
ber. The dependence of the energy resolution of a scintil- 
lation detector vs energy can be written as a = C x \AE. 
The parameter C was found to be 0.063 MeV 1/2 and 0.068 
MeV ' for the first and second detectors, respectively. 
The values of a determined from the background spectrum 
are in good agreement with measurements performed with 
a 207 Bi calibration source. The expected dispersion of the 
5.5 MeV peak due to axion absorption is a — 148 keV and 
160 keV for BGOl and BG02, respectively. 

Figure 0] shows the energy range of (4.5 — 6.5) MeV, 
in which the axion peak was expected. Since the energy 
resolutions of the detectors differ by 10%, the BGOl and 
BG02 spectra were analyzed separately. The spectrum 
measured in the range of 4.7 — 6.4 MeV was fitted by a 
sum of linear and Gaussian functions; the Gaussian peak 
position and dispersion corresponded to the desired-peak 
parameters: 



N(E) =a + bE 



-exp[ 



(go - Ef 

2CT 2 



where Eq — 5.49 MeV is the axion peak position, a = 
0.148 (0.160) MeV is the Gaussian peak dispersion, S is 
the peak area, and a and b are the parameters of the func- 
tion describing the continuous background. 

The position peak and dispersion were fixed and three 
parameters were varied, two of which described the contin- 
uous background while the third described the area of the 
desired peak. The total number of the degrees of freedom 
in the range of 4.5-6.5 MeV was 134 (139). The fit results 
for BGOl, corresponding to the minimum \ 2 ~ 139/134, 
are shown in Fig. 4. The intensity of the 5.5 MeV peak 
was found to be S = -44 ± 78 (S = -4 ± 81 for BG02). 
Combining these results, we obtain S = — 48 ± 111 for the 
axion peak area; this corresponds to the upper limit on 
the number of counts in the peak, Su m = 140 at a 90% 
confidence level |25j . 

The expected number of axioelectric absorption events 



Sabs — £N B iT<pA&Ae 



(13) 



where a Ae is the axioelectric effect cross section, given by 
expression ©; <Pa is the axion flux ([7]); Nsi = 1.88 x 
10 24 is the number of Bi atoms; T — 5.62 x 10 6 s is the 
measurement time; and e = 0.87 is the detection efficiency 
for 5.5 MeV electrons. Axion flux <Pa is proportional to the 
constant g AN , and the cross section a Ae is proportional 
to the constant g Ae , according to expressions ([7]) and 
©. As a result, the S a bs value depends on the product of 
the axion-electron and axion-nucleon coupling constants: 



9Ae 



2 

x 9 an- 



The experimentally found condition S a bs < Su m im- 
poses some constraints on the range of possible gAe x 9an 
and vnA values. The range of excluded \gAe x 5a at I values 
is shown in Fig. 5, at rriA — > the limit is \gAe x 9\n\ — 
3.2 x 10~ 9 . The dependence of \gAe x 9an\ on mA ^ s re ~ 
lated only to the kinematic factor in formulas © and (JSJ) . 
These constraints are completely mo del- independent and 
valid for any pseudoscalar particle. 

Within the hadronic axion model, gAN and rriA quan- 
tities are related by expression ([S]), which can be used to 
obtain a constraint on the gAe constant, depending on the 
axion mass (Fig. 5). For mj = 1 MeV, this constraint 
corresponds to gAe < 1-8 x 10~ 7 . Figure 5 also shows the 
constraints on the constant gAe that were obtained in the 
Borexino experiment for 478-keV solar axions [25] and in 
the Texono reactor experiment [27] , where an attempt was 
made to observe the Compton conversion of 2.2-MeV ax- 
ions produced in the n + p — > d + A reaction. 

The experiment's sensitivity to the constant gAe de- 
pends on the target mass M, specific background level B, 
detector resolution a, and measurement time T: 



9Ae/S( 9 Ae) = (MT/aB)^ 4 . 



(14) 



(12) 



The record level of sensitivity gAe ~ 10~ 8 can be achieved 
with the detector mass enlarged by an order of magnitude 
and the background level reduced by two orders of mag- 
nitude (the latter can be done by placing the setup in an 
underground laboratory). 
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Fig. 5. Range of excluded gAe and \gAe x SajvI values (above 
the curves), depending on mA, as compared to the results of the 
(1) Borexino [21] and (2) Texono [27J experiments. The allowed 
gAe and \gAe x SajvI values lie below the corresponding curves. 



6 CONCLUSIONS 

A search for the reaction of the axioelectric absorption of 
5.5 MeV axions produced in the p + d — > 3 He + 7 (5.49 
MeV) reaction was conducted using two BGO detectors 
with a total mass of 1.2 kg, located in a low-background 
setup equipped with passive and active shieldings. As a 
result, the constraints on the axion-electron coupling con- 
stant gAe < (1.8 — 9.0) x 10~ 7 for axions with masses 
0.1 < wia < 1 MeV were obtained. For axions with masses 
mA > 2m e , we calculated the solar positron flux due to the 
decay of axions into electrons and positrons:A — > er +e + . 
The measured value of the interplanetary positron flux 
allowed us to establish a new constraint on the axion- 
electron coupling constant: gAe < (1— 5) x 10~ 17 for axions 
with masses in the range of (1.2 — 5.4) MeV. 
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